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Abstract 

Recently, there has been a wide interest in the study of aggregation equations and Patlak-Keller-Segel 
(PKS) models for chemotaxis with degenerate diffusion. The focus of this paper is the unification and 
generalization of the well-posedness theory of these models. We prove local well-posedness on bounded 
domains for dimensions d > 2 and in all of space for d > 3, the uniqueness being a result previously 
not known for PKS with degenerate diffusion. We generalize the notion of criticality for PKS and show 
that subcritical problems are globally well-posed. For a fairly general class of problems, we prove the 
existence of a critical mass which sharply divides the possibility of finite time blow up and global existence. 
Moreover, we compute the critical mass for fully general problems and show that solutions with smaller 
mass exists globally. For a class of supercritical problems we prove finite time blow up is possible for 
initial data of arbitrary mass. 

1 Introduction 

Nonlocal aggregation phenomena have been studied in a wide variety of biological applications such as migra- 
tion patterns in ecological systems [16j [60] [5_TJ [29l [17] and Patlak-Keller-Segel (PKS) models of chemotaxis 
[Ml [Ml 1301 133 HO] ■ Diffusion is generally included in these models to account for the dispersal of organisms. 
Classically, linear diffusion is used, however recently, there has been a widening interest in models with de- 
generate diffusion to include over-crowding effects [501 E] ■ The parabolic-elliptic PKS is the most widely 
studied model for aggregation, where the nonlocal effects are modeled by convolution with the Newtonian 
or Bessel potential. On the other hand, in population dynamics, the nonlocal effects are generally modeled 
with smooth, fast-decaying kernels. However, all of these models are describing the same mathematical phe- 
nomenon: the competition between nonlocal aggregation and diffusion. For this reason, we are interested 
in unifying and extending the local and global well-posedeness theory of general aggregation models with 
degenerate diffusion of the form 

u t + V • (uv) = AA(u) in [0, T) x D, (la) 
v = X7K*u. (lb) 

Mathematical works most relevant this paper are those with degenerate diffusion [51 [571 EH1 (Ml EH [Ml [Ml HFJ 
and those from the classical PKS literature [Ml 123 H5J Q3] . See also [55] , 

Existence theory is complicated by the presence of degenerate diffusion and singular kernels such as the New- 
tonian potential. Bertozzi and Slepcev in [8J prove existence and uniqueness of models with general diffusion 
but restrict to non-singular kernels. Sugiyama 59* proved local existence for models with power-law diffusion 
and the Bessel potential for the kernel, but uniqueness of solutions was left open. We extend the work of 
[5] to prove the local existence of ([1]) with degenerate diffusion and singular kernels including the Bessel 
and Newtonian potentials. The existing work on uniqueness of these problems included a priori regularity 

* jacob.bedrossian@math.ucla.edu, University of California-Los Angeles, Department of Mathematics 
1 nrodriguez@math.ucla.edu, University of California- Los Angeles, Department of Mathematics 

* bertozzi@math.ucla.edu, University of California-Los Angeles, Department of Mathematics 



assumptions [3H] or the use of entropy solutions [IT] (see also [3T]). The Lagrangian method introduced 
by Loeper in [35] estimates the difference of weak solutions in the Wasserstein distance and is very useful 
for inviscid problems or problems with linear diffusion [47l [6] [24] . In the presence of nonlinear diffusion, it 
seems more natural to approach uniqueness in -ff -1 , where the diffusion is monotone (see [EI])- This is the 
approach taken in [U [5] , which we extend to handle singular kernels such as the Newtonian potential, proving 
uniqueness of weak solutions with no additional assumptions, provided the domain is bounded or d > 3. The 
main difference is the use of more refined estimates to handle the lower regularity of Vif * u, similar to the 
traditional proof of uniqueness of L 1 ni°°-vorticity solutions to the 2D Euler equations (64[|49] and a similar 
proof of the uniqueness of L 1 n L°° solutions to the Vlasov-Poisson equation [54] . 

There is a natural notion of criticality associated with this problem, which roughly corresponds to the bal- 
ance between the aggregation and diffusion. For problems with homogeneous kernels and power-law diffusion, 
K = c |x| 2 d and A(u) — u m , a simple scaling heuristic suggests that these forces are in balance if m = 2 — 2/d 
[13]. If m > 2 — 2/d then the problem is subcritical and the diffusion is dominant. On the other hand, if 
m < 2 — 2/d then the problem is supercritical and the aggregation is dominant. For the PKS with power-law 
diffusion, Sugiyama showed global existence for subcritical problems and that finite time blow up is possible 
for supercritical problems [SHJ [S5J [ST] . We extend this notion of criticality to general problems by observing 
that only the behavior of the solution at high concentrations will divide finite time blow up from global 
existence (see Definition [6]). We show global well-posedness for subcritical problems and finite time blow up 
for certain supercritical problems. 

If the problem is critical, it is well-known that in PKS there exists a critical mass, and solutions with larger 
mass can blow up in finite time [T3[M1IIS1[27J[T1[T3[T31[371[S51[TS]. For linear diffusion, the same critical 
mass has been identified for the Bessel and Newtonian potentials [15j [19] ; however for nonlinear diffusion, 
the critical mass has only been identified for the Newtonian potential [13]. In this paper we extend the free 
energy methods of |T3l [27] [TH] [14] to estimate the critical mass for a wide range of kernels and nonlinear 
diffusion, which include these known results. For a smaller class of problems, including standard PKS models, 
we show this estimate is sharp. 

The problem ([1]) is formally a gradient flow with respect to the Euclidean Wasserstein distance for the free 
energy 

F(u(t))=S(u(t))-W(u(t)), (2) 
where the entropy S(u(t)) and the interaction energy W(u(t)) are given by 

S(u(t)) = J $(u(x,t))dx, 

= \ J j u{x,i)K{x - y)u(y,t)dxdy. 

For the degenerate parabolic problems we consider, the entropy density <fr(z) is a strictly convex function 
satisfying 

$"(z) = ^M, $'(1) = 0, $(0) = 0. (3) 

See [23] for more information on these kinds of entropies. Although there is a rich theory for gradient flows 
of this general type when the kernel is regular and A-convex [50] [3] [22] the kernels we consider here are 
more singular and the notion of displacement convexity introduced in [50j no longer holds. For this reason, 
the rigorous results of the gradient flow theory are not directly applicable, however, certain aspects may be 
recovered, such as the use of steepest descent schemes [TT][T2]. Moreover, the free energy © is still the 
important dissipated quantity in the global existence and finite time blow up arguments. The free energy has 
been used by many authors for the same purpose, see for instance [57j [T5J [18j [T3] [5j [14]. For the remainder 
of the paper we only consider initial data with finite free energy, although the local existence arguments may 
hold in more generality. 
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There is a vast literature of related works on models similar to (TTJ). For literature on PKS we refer the 
reader to the review articles [331131]; see also [3T] [23 US] for parabolic-parabolic Keller-Segel systems. For 
the inviscid problem, see the recent works of [41] |4j [6j [22] . For a study of these equations with fractional 
linear diffusion see [HH33-MH]- When the diffusion is sufficiently nonlinear and the kernel is in L 1 , ^ may be 
written as a regularized interface problem, a notion studied in 55 . Critical mass behavior is also a property 
of other related critical PDE, such as the marginal unstable thin film equation [63, 7 and critical nonlinear 
Schrodinger equations [53] [37] . 

Outline of Paper. In Section £11.11 we state the relevant definitions and notation. Furthermore, we give a 
summary of the main results but reserve the proofs for subsequent sections. In Section we prove the 
uniqueness result. Local existence is proved in Section SJ3J The first result is proved on bounded domains 
in d > 2 and the second is proved on all space for d > 3. In Section 2] we prove a continuation theorem. 
The global existence results are proved in Section SJ5] Finally, in Section SJ5] we prove the finite time blow up 
results. 

1.1 Definitions and Assumptions 

We consider either D = M. d with d > 3 or D smooth, bounded and convex with d > 2, in which case we 
impose no-flux conditions 

(-VA(u) + uVIC * u) ■ v = on dD x [0, T), (4) 

where v is the outward unit normal to D. We neglect the case D = M. 2 for technicalities introduced by the 
logarithmic potential. 

We denote Dt := (0,T) x D. We also denote ||w|| p :— ||u||x,j>m) where LP is the standard Lebesgue space. 
We denote the set {u > k} := {x G D : u{x) > fc}, if S C R d then l^l denotes the Lebesgue measure and Is 
denotes the standard characteristic function. In addition, we use J fdx := J D fdx, and only indicate the 
domain of integration where it differs from D. We also denote the weak L p space by L p '°° and the associated 
quasi-norm 

H/lliP.oo = supa p A/(a) 

\q>0 

where A/ (a) = |{/ > a}| is the distribution function of /. Given an initial condition u(x,0) we denote its 
mass by J u(x, 0)dx — M . In formulas we use the notation C(p, fc, Af, ..) to denote a generic constant, which 
may be different from line to line or even term to term in the same computation. In general, these constants 
will depend on more parameters than those listed, for instance those associated with the problem such as 
K, and the dimension but these dependencies are suppressed. We use the notation / < Pl k,... g to denote 
/ < C(j>, fc, ..)g where again, dependencies that are not relevant are suppressed. 

We now make reasonable assumptions on the kernel which include important cases of interest, such as when 
JC is the fundamental solution of an elliptic PDE. To this end we state the following definition. 

Definition 1 (Admissible Kernel). We say a kernel K is admissible if K G W| ' and the following holds: 

(R) KeC 3 \ {0}. 

(KN) K, is radially symmetric, K,{x) = k{\x\) and fc(|x|) is non-increasing. 

(MN) k"(r) and k'(r)/r are monotone on r G (0, S) for some 5 > 0. 

(BD) \D 3 IC(x)\ < \x\- d -\ 

This definition ensures that the kernels we consider are radially symmetric, non-repulsive, reasonably well- 
behaved at the origin, and have second derivatives which define bounded distributions on LP for 1 < p < oo 
(see Section £| 1 . 3|) . These conditions imply that if /C is singular, the singularity is restricted to the origin. 
Note also, that the Newtonian and Bessel potentials are both admissible for all dimensions d > 2; hence, the 
PKS and related models are included in our analysis. 
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We now make precise what kind of nonlinear diffusion we are considering. 

Definition 2 (Admissible Diffusion Functions). We say that the function A(u) is an admissible diffusion 
function if: 

(Dl) A G C x ([0, oo)) with A'(z) > for z G (0,oo). 
(D2) A'(z) > c for z > z c for some c, z c > 0. 
(D3) ^ A' {z)z- 1 dz <oo. 

This definition includes power-law diffusion A(u) = u m for m > 1. Note that (D3) requires the diffusion 
to be degenerate at u = 0, however it is permitted to behave linearly at infinity. Furthermore, on bounded 
domains condition (D3) can be relaxed without any significant modification to the methods. Following [5], 
the notions of weak solution are defined separately for bounded and unbounded domains. 

Definition 3 (Weak Solutions on Bounded Domains). Let A(u) and K be admissible, and uq{x) G L°°(D) 
be non-negative. A non-negative function u : [0,T] x D — > [0, oo) is a weak solution to JTJ if u G L°°(Dt), 
A{u) G L 2 (0,T,H 1 (D)), u t G L 2 (0,T,H~ 1 (D)) and 

J J ' u<t> t dxdt = J ' u o (x)(j){0 1 x)dx + J J ' (V A{u) - uVJC * u) ■ V0 dxdt, (5) 

for all <j) G C°°(nP) such that cf>{T) = 0. 

It follows that uV/C *«6 L 2 (Dt); therefore, definition [3] is equivalent to the following, 

{u t {t),4>)= I (-VA(u) + uVlC*u) -V<j)dx, (6) 



for all test functions <fi G H 1 for almost all t G [0, T\, Above (•, •) denotes the standard dual pairing between 
H 1 and if -1 . Similarly for M. d we define the following notion of weak solution as in |S]. 

Definition 4 (Weak Solution in R d , d > 3). Let A and JC be admissible, and u G L°°(K d ) n L 1 (M rf ) be non- 
negative. Afunctionu : [0,T]xR d -)• [0,oo) is a weak solution of © ifu G L°°((0, T)x]R ci )nL oo (0, T, X 1 (R d )), 
G ^(OjT,^ 1 ^)), uVX*Me L 2 (D T ), u t G L 2 (0, r,^" _1 (]R d )), and for all test functions G H 1 ^) 
for a.e f G [0, T] © holds. 

We show below (Theorem [3]) that weak solutions satisfying Definition [3] or 0] are in fact unique. Moreover, 
we show the unique weak solution satisfies the energy dissipation inequality (Proposition [1} , 

F{u{t)) + J ' i \A'(u)Vu - mV/C * u\ 2 dxdt < T(u Q {x)). (7) 

As is important for the global theory, one could also refer to these solutions as free energy solutions, as is 
done in [13]. Uniqueness implies that there is no distinction between free energy solutions in [13] and weak 
solutions. 

Since (JlJ conserves mass, the natural notion of criticality is with respect to the usual mass invariant scaling 
u\(x) = X d u(Xx). A simple heuristic for understanding how this scaling plays a role in the global existence is 
seen by examining the case of power-law diffusion and homogeneous kernel, A(u) — u m and /C(x) = \x\ d ^ p . 
Under this mass invariant scaling the free energy ([2]) becomes, 

F(u x ) = \ dm ~ d S{u) - \ d / p W(u). 

As A — > oo, the entropy and the interaction energy are comparable if m = (p + 1) /p. We should expect 
global existence if to > (p + as the diffusion will dominate as u grows, and possibly finite time blow 

up if to < [p + l)/p as the aggregation will instead be increasingly dominant. We consider inhomogeneous 
kernels and general diffusion, however for the problem of global existence, only the behavior as u — > oo 
will be important, in contrast to the problem of local existence. Noting that |x| d ^ p is, in some sense, the 
representative singular kernel in L p '°° leads to the following definition. This critical exponent also appears 
indirectly in [46] . 



4 



Definition 5 (Critical Exponent). Let d > 3 and K, be admissible such that JC G L p ^ for some d/(d— 2) < 
p < oo. Then the critical exponent associated to K, is given by 

v + 1 

Km* = < 2 - 2/d. 

If D 2 K.{x) = 0(\x\~ 2 ) as x -> 0, then we take to* = 1. 

Remark 1. The case m* = 1 implies at worst a logarithmic singularity as x — > and if d = 2 then all 
admissible kernels have to* = 1. 

Now we define the notion of criticality. It is easier to define this notion in terms of the quantity A'(z), as 
opposed to using $(z) directly. 

Definition 6 (Criticality). We say that the problem is subcritical if 

r • f A '& 
nmint — r = oo, 



critical if 



and supercritical if 



A'(z) 

< liminf — f-4- < oo, 



liminf AM- = . 



Notice that in the case of power-law diffusion, A(u) — u m , subcritical, critical and supercritical respectively 
correspond to to > m*,m = m* and m < to*. Moreover, in the case of the Newtonian or Bessel potential, 
to* = 2 — 2/d and the critical diffusion exponent of the PKS models discussed in [58l [53 [13] is recovered. 

1.2 Summary of Results 

The proof of local existence follows the work of Bertozzi and Slepcev [8], where ([T]) is approximated by a 
family of uniformly parabolic problems. The primary new difficulty, due to the singularity of the kernel, is 
obtaining uniform a priori L°° bounds, which is overcome here using the Alikakos iteration [2]. Solutions are 
first constructed on bounded domains. 

Theorem 1 (Local Existence on Bounded Domains, d > 2). Let A(u) and K(x) be admissible. Let Uq(x) g 
L°°(D) be a non-negative initial condition, then (Jl]) has a weak solution u on [0,T] x D, for some T > 0. 
Additionally, u G C([0, T]; L P (D)) for p £ [l,oo). 

In dimensions d > 3 we also construct local solutions on M. d by taking the limit of solutions on bounded 
domains. 

Theorem 2 (Local Existence in R d , d > 3). Let A(u) and K{x) be admissible. Letu (x) G L°° (R d ) H L 1 (R d ) 
be a non-negative initial condition, then (fTJ) has a weak solution u on My, for some T > 0. Additionally, 
u 6 C([0, T}; L p (M. d )) for all 1 < p < oo and the mass is conserved. 

As previously mentioned, the free energy is a dissipated quantity for weak solutions and is a key tool for the 
global theory. 

Proposition 1 (Energy Dissipation). Weak solutions to (TTJ) satisfy the energy dissipation inequality (0 for 
almost all t > 0. 

As in [5], uniqueness holds on bounded, convex domains in d > 2 or on R d for d > 3. The proof also holds 
for more general diffusion (e.g. fast or strongly degenerate diffusion) or no diffusion at all. 

Theorem 3 (Uniqueness). Let D C M. d for d > 2 be bounded and convex, then weak solutions to (fTJ are 
unique. The conclusion also holds on R d for d > 3. 
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We also prove the following continuation theorem, which generalizes similar theorems used in for instance 
[TSJ The proof follows the well known approach of first bounding intermediate L p norms and using 
Alikakos iteration [2] to conclude the solution is bounded in L°° (for instance, see [38l l59l H3l l25l l34l [To] ) . 

Theorem 4 (Continuation). The weak solution to ([I]) has a maximal time interval of existence T* and either 
T* = oo or T* < oo and 

lim limsup \\(u — fc)+|| a-m > 0. (8) 

_ffere to is such that 1 < to < m* and lim inf 2 _>oc A'(z)z 1 ~ m > 0. In particular, for all p > (2 — to)/(2 — 
m*), 

lim llulL = oo. 

Remark 2. Note that the order of the limits in Theorem 0] is important. In fact, if the ordered is reversed 
the limit is always zero. 

For the case to* = 2 — 2/d, Blanchet et al. [13] identified the critical mass for the problem with the Newtonian 
potential, K, = cy \x\ ~ 2 , and A(u) — u m * . The authors show that if M < M c then the solution exists globally 
and if M > M c then the solution may blow up in finite time. There M c is identified as 

2 xl/(2-m*) 

(to* -l)C m *CdJ 

where C m * is the best constant in the Hardy-Littlewood-Sobolev inequality given below in Lemma [4] It is 
natural to ask the same question for more general cases. In this work we generalize these results to include 
inhomogeneous kernels and general nonlinear diffusion. First, we state the generalization of the finite time 
blow up results. 




Theorem 5 (Finite Time Blow Up for Critical Problems: to* > 1). Let D either be bounded and convex 
with a smooth boundary or D — R d . Let K. and A(u) be admissible and satisfy 

(Bl ) K(x) = c \x\~ d/p + o(\x\~ d/p ) as x -t for some c> and d/(d - 2) < p < oo. 

(B2) x ■ VK,{x) < -(d/p)JC(x) + C x for all x € R d , for some C x > 0. 

(B3) A'(z) = rriAz" 1 ' 1 + o(z m_1 ) as z oo for some to > I, A > 0. 

(B4) A(z) < (to - l)$(z) for all z > R, for some R > . 

Suppose the problem is critical, that is m — to*. Then the critical mass M c satisfies 

( 2A ^ 1/(2 - m<,) 



(m*-l)C m - 

and for all M > M c there exists a solution to {1} which blows up in finite time with ||tto||i = M. 

Theorem 6 (Finite Time Blow Up for Supercritical Problems). Let D be as in Theorem^ Let K, satisfy 
(Bl) and (B2) in Theorem\B[ and A(u) satisfy (B3) and (B4) in Theorem^ with 1 < to < to*. Then for 
all M > there exists a solution which blows up in finite time with ||uo||i = M. 

The Newtonian and Bessel potentials both satisfy these conditions with C\ = (Lemma 2.2, [S7]), and so 
the results apply to PKS with degenerate diffusion. Due to the decay of admissible kernels (Definition [1} 
condition (B2) should only impose a significant restriction on the behavior of K, at the origin. Power-law 
diffusion satisfies conditions (B3) and (B4); however, (B4) is also restrictive, for example, A(u) — u m — u 
for u large does not satisfy the condition. 

The accompanying global existence theorem is significantly more inclusive than the blow up theorems, both 
in the kinds of kernels and nonlinear diffusion considered. As in Theorem [5l the estimate of the critical mass 
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only depends on the leading order term of an asymptotic expansion of the kernel at the origin and the growth 
of the entropy at infinity. The approach used here and in [131 115] relies on using the energy dissipation 
inequality ([7]) and the continuation theorem (Theorem [4]). The third key component is an inequality which 
relates the interaction energy W(u) to the entropy S(u). For to* > 1 this is the Hardy-Littlewood-Sobolev 
inequality given in Lemma 2] In this case, the estimate of the critical mass is given by 

Theorem 7 (Global Well-Posedness for to* > 1). Suppose to* > 1. Then we have the following: 

(i) If the problem is subcritical, then the solution exists globally (i.e. = oo ) and is uniformly bounded 
in the sense u <E L°°((0, oo) x D). 

(ii) If the problem is critical then there exists a critical mass M c > such that if ||wo||i = M < M C! then 
the solution exists globally and is uniformly bounded in the sense u G L°°((0, oo) x D). The critical 
mass is estimated below in ([9]). 

Proposition 2 (Critical Mass For to* > 1). If K = c \x\ d ^ p + o(\x\ d ^ p ) as x — >• for some c > and p, 
d/[d — 2) < p < oo, then M c satisfies, 

li m £M _ 9^L cM 2-ra* = (g) 
z^oc z m 2 

If c = or liniz-^oo <$>(z)z~ m = oo then we define M c = oo. 

Remark 3. By LemmaHH tf £ G L p a T then 3 5 > C > such that V:E ' \ x \ < S > /C ( a; ) ^ C \x\' d/p - Then, if the 
kernel does not admit an asymptotic expansion as in Proposition [2l the critical mass M c can be estimated 

by, 

li m £M _ 9niL C M 2 c - m " = 0. 

z^oo z m 2 

Remark 4. Note, lim^oo $(z)z _m is always well-defined but is not necessarily finite unless 

lim sup A'(z)z 1 ^ m * < oo. 

z— >oo 

If the problem is critical then necessarily liniz-i.oo <f>(z)z _m > so there always exists a positive mass which 
satisfies ©• Moreover, if the problem is subcritical then necessarily lim^oo $(z)z~ m = oo. 

The case to* = 1 is analogous to the classical PKS problem in 2D, where linear diffusion is critical. For the 
2D PKS, the critical mass is given by M c = 8ir for both the Newtonian and Bessel potentials [HI [19]. In this 
work we treat the to* = 1 case for d > 2 on bounded domains, recovering the critical mass of the classical 
PKS, although (D3) technically requires the diffusion to be nonlinear and degenerate. The case d > 3 and 
to* = 1 is approached in [35] . but the optimal critical mass is not identified. Our estimate is given below in 
(|10[) . As above, the critical mass only depends on the asymptotic expansion of the kernel at the origin and 
the growth of the entropy at infinity. We first state the analogue of Theorem [5] 

Theorem 8 (Finite Time Blow Up for Critical Problems m* = 1). Let D be a smooth, bounded and convex 
domain and d > 2. Suppose K, satisfies 

( CI ) K,{x) = — cln |x| + o(ln |x|) as x — ¥ for some c > . 
(C2) x • V/C(x) < -c + C \x\ for all x G R d , for some C > . 
(C3) A(z) < Az for some A>0. 
Then the critical mass M c satisfies 

M c = , 

c 

and for all M > M c there exists a solution which blows up in finite time with ||wo||i = M . 

The corresponding global existence theorem includes more general kernels and nonlinear diffusion. The proof 
is similar to Theorem except that the logarithmic Hardy-Littlewood-Sobolev inequality (Lemma [5]) is used 
in place of the Hardy-Littlewood-Sobolev inequality. 
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Theorem 9 (Global Well-Posedness for m* = 1 on Bounded Domains). Suppose m* — 1 and d > 2, let D 
be bounded, smooth and convex. Then we have the following: 

(i) If the problem is subcritical, then the solution exists globally and is uniformly bounded in the sense 
u G L°°((0,oo) x D). 

(ii) If the problem is critical then there exists a critical mass, M c > 0, such that if ||uo||i — M < M c , then 
the solution exists globally and is uniformly bounded in the sense u G L°°((0, oo) x D). The critical 
mass is estimated below in (1101). 



Proposition 3 (Critical Mass for m* — 1 on Bounded Domains). If /C(x) = — cln|x| + o(ln|x|) asuO 
for some c > 0, then M c satisfies, 

Urn - ±M C = 0. (10) 

z->oo zmz Id 



If c = or lim^^oo $>(z)(zln z) 1 = oo then we define M c 



oo. 



Remark 5. By (BD) and (MN), 36, C > such that Vx, |x| < S, /C(x) < -Clnx. Therefore, if the kernel 
does not have the asymptotic expansion required in Proposition^ then the critical mass M c may be estimated 
as, 

2-s-oo zmz 2d 

Remark 6. These theorems include many known global existence and finite time blow up results in the 
literature including [551 EZ1 EFJ HI EH I3FJ HFJ ■ Our main contributions to the existing theory is the unification 
of these results and the estimate of the critical mass for inhomogeneous kernels and general nonlinear diffusion. 
In the case of the Newtonian potential Blanchet et al. showed in [T3] that solutions at the critical mass also 
exist globally. See [37JHU1H1] for the corresponding result for classical 2D PKS. 

1.3 Properties of Admissible Kernels 

Definition [T] implies a number of useful characteristics which we state here and reserve the proofs for the 
Appendix l8.3l First, we have that every admissible kernel is at least as integrable as the Newtonian potential. 



Lemma 1. Let K be admissible. Then V/C G £d/(<*-i)>°°. Ifd>3, then JC G L d /( d ~ 2 ^°°. 

In general, the second derivatives of admissible kernels are not locally integrable, but we may still properly 
define D 2 K, * u as a linear operator which involves a Cauchy principal value integral. By the Calderon- 
Zygmund inequality (see e.g. [Theorem 2.2 |56]]) we can conclude that this distribution is bounded on L p for 
1 < p < oo. The inequality also provides an estimate of the operator norms, which is of crucial importance 
to the proof of uniqueness. 

Lemma 2. Let /C be admissible and v — VIC * u. Then Vp, 1 < p < oo, 3 C(p) such that || Vw|| p < C(p)||u|| p 
and C(p) < p for 2 < p < oo. 

One can further connect the integrability of the kernel with the integrability of the derivatives at the origin, 
which provides a natural extension of Lemma [2] through the Young's inequality for L p, °°. 

Lemma 3. Let d > 3 and K, be admissible. Suppose 7 is such that 1 < 7 < d/2. Then K G L dl f h 2) '°° 
if and only if D 2 fC G L 1 ^ . The same holds for V/C G L d ^ oc 1 - >: °°. Ln particular, m* = 1 + I/7 — 2/d for 
some 1 < 7 < d/2 if and only if D 2 K, G L]£° . Moreover, m* = 1 if and only if D 2 K G L d {^'°° . 

The following lemma clarifies the connection between the critical exponent and the interaction energy. 

Lemma 4. Consider the Hardy- Littlewood-Sobolev type inequality, for all f G LP , g G L q and K. G for 
1 < p, q, t < 00 satisfying l/p+ 1/q + 1/t = 2, 



f{x)g{y)K,{x - y)dxdy 



< 



\p\\9U\K>\\l^. (11) 
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See 144V - I n particular, if {p + l)/p = to* > 1, then for all ueL'n L m , 

u{x)u(y) \x - yr d/p dxaly < C m * ||u||?- m * \\u\\^ . (12) 



Here C m * , depending only on p and d, is taken to be the best constant for which (|12p holds for all such u. 

Remark 7. It is not necessarily the case that C m * is easily related to the optimal constant in (JTTJ) . It is shown 
in |13] that C 2 _2/d is achieved for a fairly explicit family of extremals, but to our knowledge, extremals of 
([T21 have not been constructed for other values of m*. 



If to* = 1 then we will need the logarithmic Hardy-Littlewood-Sobolev inequality, as in for instance [27l [14] . 

Lemma 5 (Logarithmic Hardy-Littlewood-Sobolev inequality [20.). Let d > 2 and < / 6 L 1 be such that 
| J / In fdx\ < oo. Then, 



f(x)f(y)ln\x-y\dxdy<Mk I f]nfdx + C(\\f\\i)- (13) 
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2 Uniqueness 

We now prove the uniqueness of weak solutions stated in Theorem [3] 

Proof. (Theorem[3]) The proof follows [31 IB] and estimates the difference of weak solutions in H~ x , motivated 
by the fact that the nonlinear diffusion is monotone in this norm [5T] . To this end, if the domain is bounded, 
we define <fi(t) as the zero mean strong solution of 

A<j>(t) = u(t) - v(t) in D (14) 
V<f>(t) -v = 0, on dD, (15) 

where v is the outward unit normal of D. If the domain is R d for d > 3, we let <p(i) = —N* (u — v) where J\f 
is the Newtonian potential. In either case, by the integrability and boundedness of weak solutions u(t) and 
v(t) we can conclude (f>(t) e L°°(D T ) n C([0, T];H X ), V<j>(t) £ L°°(D T ) n L 2 (D T ) and <p t solves, 

A(j> t = dfU — d t v. 

Then since \\u(t) — = ||V0(i)||2, we will show that ||V</>(t)||2 = 0. During the course of the proof, 

we integrate by parts on a variety of quantities. If the domain is bounded, then the boundary terms will 
vanish due to the no-flux conditions (|4"1) . (jT5j) . In R d , the computations are justified as V/C * u, VA(ji), V/C * 
v, VA(v), Vcj) € L 2 {D T ). 

By the regularity of <fi(t) and the no-flux boundary conditions (1151) . Q we have possibly up to a set of measure 



zero, 



Therefore, using <f>(t) in the definition of weak solution and (|15l) we have, 



ld_ 

Ydt 



\V(j)(t)\ 2 dx = I (VA{u{t)) - VA(v(t))) ■ y<p{t)dx 



[u - u)(V/C * u) ■ \7cj)dx - J u(V£ * (u - v)) ■ V<t>dx. 
:=h+h + h. 

We drop the time dependence for notational simplicity. Since A is increasing, we have the desired monotonicity 
of the diffusion, 



h = - I (A(u) - A(v)) (u - v)dx < 0. 
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We now concentrate on bounding the advection terms. 
We follow |S]. By integration by parts we have, 

I2 = ^ J dj(j){dijK * u)dj4>dx + ^ J di4>{djK * u)dij(j)dx. (16) 



If the domain is bounded, we may apply integration by parts, 



/ di4>{djK * u)dij<j)dx = — / djjcfrdjK * udj(f>dx — / di<f){djjK, * u)di4>dx 
i,j i,3 i,j 

+ ~S~] / \di(j)\ 2 dj)C * ui/jdS, 

where 1^ is the unit outward normal to 13. As in [8], we have V/C * m ■ y < on 913 since 13 is convex and K, is 
radially decreasing, so that term is non-positive. If the domain were R d , such boundary terms would vanish. 
Therefore by integration by parts again we have, 

J di4>(djK. * u)dij(j)dx < — ~ / (A/C * it) |V0| 2 da;, 
which together with (fl6|) implies, 



^2 < y |£» 2 /C*u| V0| 2 dx. 
By Holder's inequality, Lemma [2] and V</> 6 L°°(Dt) for p > 2, 

O-l)/}) 



y |Z) 2 /C * w| |V0| 2 da- < \\D 2 K,*u\\ p (j \V4>\ 2p/{p ~ 1] dx 



/ r- \ (p-i)/p 

^p[i«iip[[v0ii^ (j iv0i 2 dxj 



J\V4>\ 2 dx\ , (17) 

where the implicit constant depends only on the uniformly controlled L p norms of u and w. 

As for I3, we compute as in [8]. By the computations in the proof of Lemma [2] we may justify integration by 
parts on the inside of the convolution, that is, 

||£ J d i K{x-y)dj i ff>dxh < l|V0|| 2 . 
3 

which by Cauchy-Schwarz implies, 

/ 3 <IMU|V0|| 2 . (is) 

Letting 77(f) = J |V0(f)| 2 dx, (fT7|) and (|18l) imply the differential inequality, 

where C again depends only on the uniformly controlled LP norms of u, v. The differential equality does not 
have a unique solution, but all of the solutions are absolutely continuous integral solutions bounded above by 
the maximal solution 77(f). By continuity, for t < 1/C the maximal solution is given by 77(f) = (Ct) p , hence, 

77(f) < 77(f) - (cty. 
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For t < 1/(2(7) we then have 

m < fj(t) < 2-*, 

and we take p — ► oo to deduce that for i 6 [0,1/(2(7)), rj(t) — 0, therefore the solution is unique. This 
procedure may be iterated to prove uniqueness over the entire interval of existence since the time interval 
only depends on uniformly controlled norms. □ 



3 Local Existence 

3.1 Local Existence in Bounded Domains 

Let A(z) be a smooth function on R + such that A (z) > r\ for some 77 > 0. In addition, let v be a 
given smooth velocity field with bounded divergence. Classical theory gives a global smooth solution to the 
uniformly parabolic equation 

u t = AA(u) - V • (uv) (19) 
(see [IS]). The solutions obey the global L°° bound 

||«ll L oc (D) <||«o|| i o. (D) e l|(V ' i °-^-^ t . (20) 

We take advantage of this theory to prove existence of weak solutions to |T]) by regularizing the degenerate 
diffusion and the kernel. Consider the modified aggregation equation 

u e t = AA e «) - V • (u e {VJ e 1C * O) , (21) 

with corresponding no-flux boundary conditions ((4]). We define 

A e {z) = { a' e (z)dz, (22) 



where a' e (z) is a smooth function, such that A'(z) + e < a' e (z) < A'{z) + 2e, and the standard mollifier is 
denoted J t v. We first prove existence of solutions to the regularized equation ([2"T]) , this is stated formally in 
the following proposition. 

Proposition 4 (Local Existence for the Regularized Aggregation Diffusion Equation). Let e > be fixed 
and uq(x) € C°°{D), then (|21[) has a classical solution u on Dt for all T > 0. 

We obtain the proof of Proposition |4] directly from Theorem 12 in [8]. The proof requires a bound on 
||VA c ||i2 ( - I ) T ), for some T > 0. We state this lemma for completeness but reference the reader to [8] for a 
proof. 



Lemma 6 (Uniform Bound on Gradient of A(u)). Let e > be fixed and u e 6 L°°(Dt) be a solution to ([2Tj) . 
There exists a constant C = C(T, ||Vj7e/C * u e \\ Loa ^ , Hu 6 ^) such that: 

l|VA £ K)|| i2(I)T) <C. (23) 
Remark 8. The estimate given by (f23| is independent of e. 

Proposition^] gives a family of solutions {it e } e >o- To prove local existence to the original problem ([1} we first 
need some a priori estimates which are independent of e. Mainly, we obtain an independent-in-e bound on 
the L°° norm of the solution and the velocity field. This is the main difference in the local existence theory 
from [8]. Due to the singularity of the kernels significantly more is required to obtain these a priori bounds. 
We first state a lemma, due to Kowalczyk [38] and extended to d > 2 and M. d in 18 . The proof is based on 
the Alikakos iteration. 

Lemma 7 (Iteration Lemma [381 118]). Let < T < 00 and assume that there exists a c > and u c > such 
that A' (u) > cfor all u > u c . Then if ||V/C * < (7i on [0,T] then WuW^ < C 2 {C 1 ) max{l, M, |KIL} on 
the same time interval. 
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Lemma 8 (L°° Bound of Solution). Let {u e } e>0 be the classical solutions to (l2~Tj) on Dt, with smooth, non- 
negative, and bounded initial data J e UQ. Then there exists C — C(||ito||i, ||ito||oo) cmd T = T(||ito||i, ||uo||p) 
for any p > d such that for all e > 0, 



\\u<(t)\ 



L°°(D) 



< c 



(24) 



for all t G [0, T] 



Proof. For simplicity we drop the e. The first step is to obtain an interval for which the LP norm of u is 
bounded. Following the work of |34j we define the function ui = (u e — fc)+, for k > 0. Due to conservation 
of mass the following inequality provides a bound for the LP norm of u given a bound on the LP norm of u k , 



Ni^c^XF^Hi + HliP- 

We look at the time evolution of ||wfc|L an d make use of the parabolic regularization (|22p. 
Step 1: 

-t, =P «£ -1 V • (VA £ (u) - uVJ t K * u) dx 



(25) 



-fcH / A' (it) 



Vu 



P/2 



da; + p{p — 1) / u? 1 Vufc • V JJZ * u dx 



kp { p^)JuVVu k .VJ^u d x, 



where we used the fact that for I > 

u(u k y = {u k y +i + ku 

Hence, integrating by parts once more gives 



(26) 



|kii;<^M'(i<' 



p/2 



dx — (p — 1) J u p ,AJ' e )C*udx — kp J u p k 1 AJ t tC*udx 



< -C(p)jA'(u) 

< -C(p)[a'(u) 



Vu 



p/2 



dx + C(p) \\u k \\ p p+1 \\AJ e K * u\\ p+1 + C(p)k Wukf- 1 \\AJ e K * u\\ p 



Vu 



p/2 



dx 



C{p)(\ 



iik 



ip+i 
lp+i 



»||ffi)+t7(p)fc(||u fc ||J + ||«||5 



In the last inequality we use Lemma [5] Now, using (l23t we obtain that 



-\\u k \\ p p dx<-C(p)lA'(u) 



Vu 



P/2 



da: + C(p) + C(p, k) \\u k f p + C(p, k, M). 



An application of the Gagliardo-Nirenberg-Sobolev inequality gives that for any p such that d < 2(p+ 1) (see 
Lemma flfl in the Appendix): 



,P/2 



Oil 



where a% = d/p, a 2 = 2(p+ 1) — d. From the inequality a r b^ 1 r ^ < ra+ (1 — r)b (using that a = S ||u p ^ 2 || ^1,2 
and r = ct\/2) we obtain 



\\<tlZ^\\< 



ft , r(5 2 



•r<i a ||u||5. 
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Above /3i, /?2 > 1- For k large enough we have that A'{u) > c > over {u > k}; hence, if we choose S small 
enough we obtain the final differential inequality: 

j t 1|< < C(p) \\u k f p > + C(p, k, rS) \\u k \f p + C(p, k, Huollx). (27) 

The inequality (|2"T)) in turns gives a T p — T(p) > such that ||tife|| is bounded on [0,T p ]. Inequality (l23t 
gives that ||u|| p remains bounded on the same time interval. Next we prove that the velocity field is bounded 
in L°°(D) on some time interval [0, T\. This then allows us invoke Lemma [7] and obtain the desired bound. 

Step 2: 

Since V/C G Lj oc and V/Cl R d\ Bl ( ) G L q for all q > d/(d — 1) (by Lemma [J), we have for all p > d/(d — 1), 

\\v\\ p = || V/C * u|| p < ||V/Cl Bl (o)||i||«||p + ||V/Cl Rrf \ fll(0) || p M. 

By Lemma Owe also have, for all p, 1 < p < oo, 

||W|| P = \\D 2 lC*u\\ p < \\u\\ p . 

By Morrey's inequality we have v G L°°(Dt) by choosing some p > d and invoking step one, and Lemma [7] 
concludes the proof. Note that the bound depends on the geometry of the domain through the constant on 
the Gagliardo-Nirenberg-Sobolev inequality (Lemma [17]). However, this constant is related to the regularity 
of the domain, and not directly to the diameter of the domain. □ 

In addition to the a priori estimates the proof of Theorem [JJ requires precompactness of {tt e } e >o in L^^Dt)- 

Lemma 9 (Precompactness in L 1 (ri-r))- The sequence of solutions obtained via Proposition^ {u c } e> o, 
which exist on [0, T], is precompact in L 1 {Dt)- 

The proof of Lemma [9] follows exactly the work in [8] . The key is to prove that the sequence satisfies the 
Riesz-Frechet-Kolmogorov Criterion. This relies on the fact that \\A(u e )\\i J 2(o,T:H 1 (D)) < C uniformly. 

Proof. (Theorem [lj For a given e > 0, if we regularize the initial condition uf)(x) = J<:Uq{x), Proposition |4] 
gives a solution u e to (|21[) . Furthermore, the proof of Proposition 2] and Lemma [5] gave uniform-in-e bounds 
on m e («)|| L 2( t H 1 (£>))' II u£ IIl=°(,d t )> ano - W u t \\l 2 (o t h- 1 (d))- By Lemma[8j all solutions exist on [0, T], with 
T independent of e. Also, recalling that A e (z) > A(z) and a' e (z) > A'(z) gives that 

m( ue )lli2(o,r,ii' 1 (o)) ^ ^ 

where C is independent of e. Since L 2 (0,T, H 1 (D)) is weakly compact there exists a p such that some 
subsequence of {u e } f>0 converges weakly, i.e A(u ej ) — p in L 2 (0, T, H 1 (D)). Precompactness in L 1 implies 
strong convergence of u ej to some u G L 1 (I?r); therefore, A(u) = p. In fact, the L°°(Dt) bound on u ej gives 
strong convergence in L p (Dt), for 1 < p < oo, via interpolation. Also, Young's inequality gives 



||^V^/C * - uV/C * u\\ lHDt) < ||u|| iDO(lM || VJ £j /C * - V/C * «|| l1(JJt) 

+ ||VJ £j X*^|| LOC(Dr) ||u e * - ulUi^) 

< (|I«IIloo (0t) llv/c|| LL; + ||vj e3 /c * ^ |U (1?r) ) IK J - u\\ l1{Dt) 

+ NUK' HoollV^/C - V£|| £ i oe . (28) 

Therefore, by interpolation u satisfies ([5]). Furthermore, we obtain that u G C([0, T]; H^ 1 (D)). To prove 
that u(t) is continuous with respect to the weak L 2 topology one uses standard density arguments. Since 
D is a bounded, u is therefore also continuous in the weak L 1 topology. To prove continuity in the strong 
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L? topology we define F(z) = J Q Z A(s)ds and show that it is continuous in the strong L 1 topology. Indeed, 
Lemma Q3] in the Appendix, see [5] for a proof, gives 



lim 

h->o 



(F(u(t)) - F(u(t + h))) dx 



t+h 



= lim / < u t ,A(t) > dr. (29) 



Recall that \\A(u)\\ l ^ (Dt) < A(||w|| L ^ (I>r) ) and so A{u) G L 2 (0, T, H~ 1 (D)). Therefore, the left hand side of 
(|29|) goes to as h 0. Now, we can invoke Lemma H"5l in Appendix, [8 , to obtain that u G C([0, T]; L 2 (D)). 
Using interpolation the L°° bound of u gives that u G C([0, T]; L P (D)), for 1 < p < oo. □ 

3.2 Local Existence in M d 

Now we consider solutions to ([1]) in R d for d > 3. We obtain such solution by taking the limit of the solutions 
in balls centered on the origin with increasing radius n, denoted by B n . 

Proof. (Theorem [2]) Let B n be defined as above and consider the truncation of the initial condition on B n , 
i.e. Uq = 1b„uq- By TheoremQ] we have a family of solutions {it n } n >o on B n for all t G [0, T]. Define a new 
sequence, {u n } n >0i where u n is the zero extension of u n . The previous work for bounded domains gives the 
uniform bounds 

\\un\\ L ^ (K) <Cx, (30) 
\\VA{u n )\\ L2[K) <C 2 . (31) 

The bounds may be taken independent of n since the constant in the Gagliardo-Nirenberg-Sobolev inequal- 
ity, Lemma [171 does not depend directly on the diameter of the domain and may be taken uniform in n — > oo. 

Therefore, there exist u,w G i 2 (R^) for which u n — u and VA(u n ) — ^ w in i 2 (My). Furthermore, (1501) 
implies M| £ oo( R d j < C\. Precompactness of {u%} t>0 in L 1 (_B„) for fixed n > and Theorem 2.33 in [I] gives 

that {u n } n >o is precompact in Lj oc (M.^). Therefore, up to a subsequence, not renamed, u n — > u in L i 1 oc (R^); 
thus, u; = VA(«). Also, the L°° bound gives that u n — > u in L^ oc (IR ti ) for 1 < p < oo. 

In addition, we have the estimate 

||G„V/C * Un|| z ,2 (K i ) < || V/C * Un|li=o (H i) l|u™|| L 2 (R d ) • (32) 

Therefore, we can extract a subsequence that converges weakly to some W\ G L 2 (M.^). Since u1b„ G 
L oo (0,T,L 1 (M d )) and ul Bn / u by monotone convergence u G L°°(0, T, L 1 (Rj.)). Once again, from the 
estimates performed in the bounded domains u n X/K. * u n —> mV/C * u in Lj oc (M.^). Therefore, we can identify 
w\ = uS/KL * u. 

We now show that u G C([0, T]; L} oc (R d )), which we know to be true, implies that u G C([0, T];L 1 (D)). 
Let t n — > t G [0, T] then for all R > we have, 



\u(t n ) — u(t)\ dx = I \u(t n ) — u(t)\ dx + / \u{t n ) — u(t)\ dx. (33) 

J Br JR d \Bn 

The first term on the right hand side of (|3"3")l can be bounded by e/2, provided n is chosen large enough, 
since u G C([0, T]; L} oc (R d )) . To bound the second term we first show that A(u) G L 1 ^). By (D3) we 
can deduce lim z ^o A(z)z^ 1 = 0. Then, for k > there exists some < Ck < oo such that if z < k then 
A(z) < Cz. Hence, 



A(u)dx = I A(u)dx + / A(u)dx 

J{u<k} J{u>k} 

< CM + A(||«|| 00 )A u (fc) < oo. 
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Therefore, ||.A(u) Hl^r;*, ) < C(M, ||u|| 00 )T. Now, let w(x) be a smooth radially-symmetric cut-off function 
with w(x) = for |x| < 1/2 and w(x) = 1 for \x\ > 1. Then consider the quantity, Mjj(t) = J uw(x/ R)dx. 
Then formally, 

j t M R {t) = Juv ■ {Vw){x/R)dx + jp J A(u)(Aw)(x/R)dx. 
Estimating terms in L°° gives, 

jM R (t) < — ^- + - 2 J A{u)dx. 

Formally, then 

M R {t) < M R (0) + M\\v\\ LK{0>t) , LX) R- x + \\A(u)\\ Limt)xRd) R- 2 . (34) 

Since A e ^((O,*) x R d ) and M R (0) — > as R -¥ oo, by choosing i? sufficiently large, the last term of ([33| 
can be bounded by e/2. Hence, implies that u G C([0, T]; L 1 (R d )). Furthermore, via interpolation we obtain 
that u e C([0,T]; L p (R d )) for 1 < p < oo. 

Conservation of mass can be proved similarly using a cut-off function w(x) = 1 for \x\ < 1/2 and w(x) = 
for \x\ > 1, see the proof of Theorem 15 in [8] for a similar proof. □ 

We are left to prove the energy dissipation inequality ([7]). As expected, the approach is to regularize the 
energy and take the limit in the regularizing parameters. 



Proof. (Proposition [T| Define 



h( u )= r^-ds, 



then $(it) = Jq h{s)ds. The regularized entropy is defined similarly with a' e (u), as defined in (|22|) . taking 
the place of A'{u). Given a smooth solution u e to (|2ip one can verify, 

T € (u e {t))+ f ( '-L|^(u e )Vu £ - u e V JX * u e \ 2 dxdr = J" £ (u e (0)). (35) 

Here J- e (u(t)) denotes the free energy with the regularized entropy and kernel. Once again we take the limit 
e approaches zero to obtain ([7]). We first show that the entropy converges. 

Step 1: The parabolic regularization gives 

h(z) + elnz < h e (z) < h(z) + 2e\nz for 1 < z, 
h(z) + lelnz < h e (z) < h'(z) + elnz for z < 1. 

Therefore, writing <£>(u) = J Q h(s)ds + J*" h(s)ds one observes that 

- 2e < $ e (u) < + 2e(ulnu)+. (36) 
This will allow us to show convergence of the entropy. In fact, 

OH < 2e f (1 +M e lnw e ) + dx+ ||*||c±([o,|MU) / 



\u e — u \ dx. 



<2e(|I>| + ||lii« e || 00 K||i) + C||« e -tt|| 



Conservation of mass, boundedness of smooth solutions, and precompactness in L\ oc imply there exists a 
subsequence, such that as e,j — > 0, 



J <S> ej (iij)dx -> J $(u)dx. 
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Step 2: To show convergence of the interaction energy we need that for a.e t G (0, T) 

u e (t)J e IC*u t (t)dx -> u(t)K * u(t)dx. (37) 



Since /C € Lj oc (D) we know that ||/C * u\\ Lao is bounded; hence, replacing V/C with /C in (f^Hj) gives the desired 
result. Finally, we are left to deal with the entropy production functional. 

Step 3: From Lemma 10 in [23], 

/- |A'(w)Vu-uV/C* M | 2 da; <liminf / — \a' e {u e )\7u e - u £ VJ e )C * u e f dx. (38) 

J U e->0 J U e 



We also note that this was proved in [8]. The proof of ([38]) relics on a result due to Otto in [52], refer to 
Lemma [Tol in the Appendix. In our case, u e G L 1 (Dt) and J e = VA £ (u £ ) — u e VJC * u e G L\ oc {Dt)- Further- 
more, up to a sequence not renamed, u e u G L 2 and J e — 1 J in L 2 , therefore, we can apply Lemma [TBI 

For the energy dissipation estimate in M. d we again consider the family of solutions {u r } to (UJ on B r 
(for simplicity let it r denote the zero-extension of the solutions). Since w n (0)ls„ u(0) by monotone con- 
vergence we obtain that F{u n (Q)) — > F(u(0)). Noting that JC G L d ^ d ^ 2 ^ allows us to make a modification 
to (|32|) and obtain that u n K, * u n — 1 u/C * m in L 2 (K^). Furthermore, (f3T)) implies that it ra /C * u n — > uK. * u 
in L\ oc . We are left to verify the uniform integrability over all space. First note that Morrey's inequality 
implies 

l|K*6 n ||^<||vx:Hloo + ll£*^ll P 

< IIV/C * uH^ + \\)C\\ Ld /v-2),oc \\u n \\ dp/{d+2p) . 
Hence, taking p sufficiently large we obtain that /C * u n is bounded in L oc (Dt)- Therefore, 



U n fC * U n dx < \\JC * UnWog / u n dx. 
l \B k JR d \B k 

This fact along with ([34]) gives that for any e > there exists a fc e sufficiently large such that for all k > k e 

u n fC * u n dx < e. 

>\B k 

This gives convergence of the interaction energy. The result follows from the weak lower semi-continuity of 
the entropy production functional and J $(u)dx in L 2 . □ 

4 Continuation Theorem 

Continuation of weak solutions, Theorem is a straightforward consequence of the local existence theory 
and the following lemma, which follows substantially the recent work in [131 1381 115) . This lemma provides a 
more precise version of Lemma [5] and has a similar proof. 

Lemma 10. Let {u e } e>Q be the classical solutions to (|21[) on Dt, with non-negative initial data J e UQ. 
Suppose there exists T , < T < oo, such that 

sup lim sup \\(u e — fc) + || 2-m =0, (39) 

e>0 fc^oo te (o,T ) 2 ~ m * 

where m is such that 1 < m < m* and lim inf z ->oo A 1 (z)z 1 ^ m > 0. Then there exists C = C(M, |[uo||oo) such 
that for all e > 0, 

sup ||u e (i)||oo < C. 
te(o,T ) 

In particular, if Tq = oo, then {it e } £>0 ar £ uniformly bounded for all time, and therefore the weak solution 
u(t), is uniformly bounded for all time. 
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Proof. (Lemma llOp Let q = (2 — m)/(2 — m*) > 1. It will be convenient to define 7, 1 < 7 < d/2 such that 
m* = 1 + I/7 — 2/d. We first bound intermediate LP norms over the same interval, (0, To). Then we use 
Morrey's inequality and Lemma [7] to finish the proof. 

Step 1: 

We have two cases to consider, m* = 2 — 2/d and m* < 2 — 2/d, which occurs if D 2 K £ L/£° f° r 7 > 1 
(LemmaHJ. In the former we show that for any p £ (q, 00) we have u e (t) uniformly bounded in L°° (0, To; L p ). 
In the latter case we only show that for q < p < 7/(7 — 1) we have u e (t) uniformly bounded in L°° (0, T ; L v ). 
In either case, this is sufficient to apply Lemma [7] and conclude the proof. 

Let fc > be some constant to be determined later and let Uk = (u — fc)+. We have dropped the e and time 
dependence for notational convenience. By conservation of mass and (|25[) . it suffices to control ||ufc|| P for any 
fc > 0. Thus, using the parabolic regularization, (1221) . and (l25t we obtain 

j t \\u k \\ p p <-p(p~l) J u p k r 2 A'(u)\Vu\ 2 dx+p( P ~l) Jiu^r 1 + ku p k - 2 )Vu- J t VK*udx. 

Then, 



dt 



\\u k \\l<-A{p-l) / A'(u) 



Vu 



p/2 



dx 



((P- 



l)u^ + kpu v k )J' e AK. * udx. 



(40) 



Since the constants are not relevant, we treat the cases together only noting minor differences when they 
appear. Ifm = 2 — 2/d we may use Holder's inequality and then Lemma [2] to obtain a bound on the first 
term from the advcction: 



u| J^A/C * udx 



< 



p ,k IKIIp+ilMlp+i 



On the other hand, if 7 > 1 we have from the generalized Hardy-Littlewood-Sobolev inequality (fTTj) (Lemma 



vP h J t AA^ * udx 



< 



p,k ||«fc||Sp||«||t + C(M)||« fc ||P 



with the scaling condition 1/a + l/t + 1/7 = 2. Choosing t — ap implies that 

I2-I/7 



1 + 1/p 



(41) 



Notice that from our choice of p then 1 < l/p+ 1/7; thus, 1/a < 1. Note that in the case when m = 2 — 2/d 
then t = ap = p + 1. Thus we estimate the advection terms, 



u v h J e tJZ * udx 



< 



P,K ||«fc||Sp||«||ap + C(M)||u fc ||* 



ED < 



\Uk 



I n /> 

ip+i 



I ap 



C{M)\\u k \\l 
C(M)\\u k \\P p + C(k,M). 



(42) 



The lower order terms in the advection can be controlled using Holder's inequality and Lemma [2j 



p-i 



JnAK * udx 



m <\\u k \\l + C(k,M). 



(43) 

We now aim to compare the dissipation term in (j40|) with the estimates (|42|) and (|43|) . We use the Gagliardo- 
Nirenberg-Sobolev inequality (Lemma fT7|). 



||Wfc||ap 



< 



i^iirn4 p+m " 1)/2 n^. 



(44) 
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with 



and 



2d t (p - q/a) 

p \q(2 -d) + dp + d(m - 1) 



Q-2 = 1 - + m — l)/2 > 0. 
By the definition of g and (|4"Tj) we have that, 

ai(p + l)/2 = l, 

which implies, 



NII& 1 < \\u k 



,a 2 (p+l) 
I? 



Vu 



p/2 



da; - 



p+m— 1 



dx 



(45) 



(46) 



If d = 2 then necessarily m = m* = 1 and this inequality will be sufficient. However, for d > 3, more work 
must be done. Define, 



I = / u 



Vu 



p/2 



dx. 



Then, for /3i < a x and (p + m - l)/3i/2 < 1, 



2d(l-?/(p + m-l)) 



/?1 g(2 - d) + dp + d(m - 1) ' 
and /?2 = 1 — /?i(p + m — l)/2 > 0, we have the following by Lemma [P7l 

(p+m-l)ft/2 



p+m— 1 



dx 



< 



(p+m-l)/3i/2> 



Therefore, by weighted Young's inequality for products, 

^-'dx < \\ Uk £ p+m - 1)fj2 + + |K||y\ 



(47) 

for some 70 > 0, the exact value of which is not relevant. Putting (|46|) and (|47l) together implies, 

IM&ZnhkhV + CilHlh), (48) 

where Viz) denotes a polynomial such that V(z) — >• as z — > 0. By definition of m, 3 (5 > such that for 
fc sufficiently large then u > k implies A'(u) > 5u m ~ l . Therefore, combining (J40J) with (|l8 ]) .(|32" ]) and (l43t 
implies, 



-|M£<-C(p)<5 / u 



m— 1 



Vu 



p/2 



dx + C(p)\\u k f+ 1 



-C(M,p)\\u k \\P + C(k,M,p) 



< 



P+i 



P(IKH«) 

+ C(M,p)||« fc ||J + C(*,M,p J ||« fc y. 
By interpolation against L , conservation of mass and a > 1 we have 

KII^mI + IKII^ 1 . 

Therefore, by assumption (f3"5)) we may choose k sufficiently large such that there exists some 77 > which 
satisfies the following for all t £ (0,T ), 



d 
dt 



|u*||P<-»7||u fc ||P + C(fc,M,p,||u fc y. 
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It follows that ||ufc|| p is bounded uniformly on (0,Xb). 
Step 2: 

The control of these V norms will enable us to invoke Lemma [7] and conclude u e (t) is bounded uniformly 
in L°°{D To ). Since V/C £ L] oc and V/Cl ffi d\ Bl ( ) £ L q for all q > d/(d — 1) (by Lemma|l}, we have for any 
q > d/(d - 1.) 

\\v\\ q - \\VJC *u\\ q < ||V/Cl Bl(0 )||i||u||, + ||VO R , XBl(0) || g M. 

If 7 > 1, then we may choose q £ (d/(d — l),7/(7 — 1)], since in this case necessarily d > 3. Otherwise we 
may choose q > d/{d — 1) arbitrarily. Then, step one implies v £ L°°((0, To); L q ). If 7 > 1 then, noting that 
Definition [1] implies D 2 K,l^d\ Bl ^ G £ 9 f° r a ll 9 > 1, 

||W|U+i = \\D 2 1C * tiUd+i < \\D 2 Kl Bl (o) Hi",.- Hull,, + ||VO K ,\ Bl(Q) || d+ iM, 
for p = 7(d + 1) /(d(7 - 1) + 27 - 1). Note that 

1< P=T *L< ? 



d(7- 1) + 27- 1 ~ 7-1 

On the other hand, if m* = 2 — 2/d then the above proof shows that u e (t) is bounded uniformly in 
L°°((0,To);T p ) for all p < 00. Therefore, by Lemma [5] we have ||W|| P < ||u|[ p < 1, for all 1 < p < 00. 
In either case, this is sufficient to apply Morrey's inequality and conclude that ||w||oo is uniformly bounded 
on (0,To). By Lemma we then have that u £ is uniformly bounded in L°°(Dt ) and we have proved the 
lemma. As in Lemma |51 the uniform bounds depend on the domain but not it's diameter. □ 

Remark 9. The proof of this lemma directly implies global well-posedness in the subcritical case since (f3l?| 
is only necessary in the critical and supercritical cases. Moreover, in the critical case, one may prove directly 
that there exists some Mo such that if M < Mq the solution is global. However, Mq will generally depend 
on the constant of the Gagliardo-Nirenberg-Sobolev inequality, as in [SHI 1111 EH] • As discussed in the recent 
works of [131 [H] , the use of a continuation theorem will allow for a more accurate estimate of the critical 
mass through the use of the free energy. 

Proof. (Theorem [4]) Suppose, for contradiction, that the weak solution cannot be continued past T* < 00 
and (jSJ) fails. As the regularized problems are bounded, this implies the hypotheses of Lemma ITD1 are satisfied 
on (0,T,t), and therefore sup e>0 sup tg ( 0jTi -) ||u £ (i)llp ^ V as * /* 1* f° r some p >q and rj > 0. By the proof 
of Lemma El for any ij > there exists a r = t(t],M) > such that if ||ito|| p < t] then ||u e || p < C for all 
e > 0. Therefore, we may choose some t n < T+ such that r satisfies t n + r > T* and, by Theorems Q] and 
[2 we construct a solution u(x,t) on the time interval [t n ,t n + t). By uniqueness, u(x, t) — u{x,t) a.e. for 
t £ [t n ,T+)\ hence, it is a genuine extension of the original solution u(x,t). However, it exists on a longer 
time interval which is a contradiction. □ 

5 Global Existence 

We now prove Theorem [7] We first note that the entropy is bounded below uniformly in time, which is a 
consequence of assumption (D3) of Definition [2] 

Lemma 11. Let u(x,t) be a weak solution to Then, 

<P(u(t))dx > -CM. 



Proof. Let h(z) = f* A'(s)s~ 1 ds. By Definition dj (D3), for z < 1, 

h(z) > -C > -00. 
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Therefore, 



<f>(u)dx =11 h{z)dzdx > / l{ u <i} / h(z)dz + / h(z)dzdx. 



>- J l {u < 1} Cu - l {u > 1} Cdx 
> -2C||u||i. 

where the last line followed from Chebyshev's inequality. □ 

5.1 Theorem O m* > 1 

Proof. (Theorem [7]) We only prove the second assertion under the hypotheses of Proposition [21 as the 
subcritical case follows similarly. By the energy dissipation inequality ([7]) we have for all time < t < T+, 

S(u(t)) - W(u(t)) < J-(no) := F . (49) 

We drop the time dependence of u(t) for notational simplicity. By the assumption on /C, Ve > 0, 3 S > 
such that \JC(x)\ < (c + e) \x\~ d/p for \x\ < S. By Lemmal we have, 

J <P(u)dx-±C m <M 2 - m *(c + e)\\u\C: ^Fo + ^ms^W^M 2 , 

By © and M < M c , there exists e > small enough and a, k > such that 

$(z)z- ,n * - -C m *M 2 - m * (c + e) > a > 0, for all z > k. (50) 

By Lemma [TT] we have, 

/ u m * (<5>{u)u rnt - \c m *M 2 ~ m * (c + e]\dx- \ [ C m *M 2 - m * (c + e) u m Ux < F + C{5, M), 
J{u>k} V 2 J 2 J{ u<k j 

and by dSHJ, 

a [ u m *dx - -C m *M 2 ~ m * (c + e) / u m * dx <F + C(M, S). 

J{u>k} 2 J{u<k} 

By mass conservation we have that ||u|| m * is a priori bounded independent of time and Theorem [4] and 
Lemma [TOl implies global existence and uniform boundedness. □ 

5.2 Theorem M m* = 1 

The proof of Theorem [5] follows similarly, but requires the logarithmic Hardy-Littlewood-Sobolev inequality 
(Lemma O as opposed to Lemma 0] 

Proof. (Theorem [9]) 

We only prove the second assertion under the hypotheses of Proposition [31 as the subcritical case follows 
similarly. We will again use Theorem 2] and prove 



sup / (u\nu)+dx < oo. 

\(0,oo) J 



te(o 

By the energy dissipation inequality we again have (|4"5)) . By the assumptions of Proposition [3J for all 
e > there exists S > such that, 



$(u)dx+ (c+e)- / / u(x)u(y)\n\x - y\dxdy < C(F ,S, M). 

x — y | < S 



2 
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By D bounded, the logarithmic Hardy-Littlewood-Sobolev inequality ([13")) implies, 

J §{u)dx- (c+ e)— y ulnuefo < C(F , 5, M, dianxD). 

Choosing fc > large and recalling Lemma ITT1 implies 

f ( $(u) M\ f 

/ ulntt — — - (c+e) — dx- (c + e) / wlnuda; < C(F , 5, M, dianxD). 
./{«>£} \ umM 2 «/ J{n<fc} 

As in the proof of Theorem[7l by conservation of mass, ([TO]) and M < M c , we may choose e > small enough 
and k large enough such that 



/ ulnudx < C(Fo,M, diamD). 

J{u>k} 



□ 



6 Finite Time Blow Up 

In this section we prove Theorem [6] and Theorem We prove Theorem [6] as it is somewhat easier, though 
the technique is the same as that used to prove Theorem [5] 

6.1 Supercritical Case: Theorem [6] 

For Theorem |6] we state the following lemma, which provides insight into the nature of the supercritical cases. 
The proof and motivation follows [12] . 

Lemma 12. Define J>m = {ael'fl L m : u > 0, ||u||i = A/}. Suppose JC satisfies (Bl) and A(u) sat- 
isfies (B3) for some m > 1,A > 0. Suppose further that the problem is supercritical, that is, m < m* . 
Then 'vafy M J- = — oo. Moreover, there exists an infimizing sequence with vanishing second moments which 
converges to the Dirac delta mass in the sense of measures. 

Proof. Let < 6 < 1, a = d/p. Then by Lemma 0] there exists h* such that, 

1 1 h*(x)h*(y) \x - y\~ a dxdy 



9C m * < 



\\h*\\l- m *\\h 



<G m ,. (51) 



We may assume without loss of generality that h* > 0, since replacing h* by \h*\ will only increase the value 
of the convolution. By density, we may take h* £ and therefore with a finite second moment. 

Let (j, = ll/i'll^Af- 1 ^, A > and h x {x) = \ d h*(\fix). First note, by (B3), Ve > 0, 3R > such that, 



J <$>{h x )dx = J J X J ^-^dzdsdx 



< I I / (mA + e)z m - 2 dz+ / —^Idzdsdx 

'o Jr Ji 



z 



<^-\\hx\\2 + C{R)\\hx\\i. (52) 

TO — 1 



By (Bl) and h* e C c °°, Ve > 0, 3X > sufficiently large such that, 

p,- 2d+a \ a 



- W(t) < -(c - e)!=—i / / h*(x)h*(y) \x - y\~ a dxdy. (53) 

Combining (f53 |) . (l52j) with (|5T|) and Lemma |H we have for X,R sufficiently large, 
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\dm-d-\/f (-i / \\h*\\ \- 2 + a / d 

w < (ro A _ 1) ||^| U ^+e)re - n<>- e)?? (i>yi) n^ir- h^iis: 

+ C(R) tJ ,- d \\h*\\ 1 . 

By supercriticality, we have a — dm* — d > dm — d, and so for e < 0, we take A — > oo to conclude that for 
all values of the mass M > we have inf T = — oo. Moreover, since h* € C^°, the second moment of /ia 
goes to zero and h\ converges to the Dirac delta mass in the sense of measures. □ 

Proof. (Theorem [6]) We may justify the formal computations for weak solutions using the regularized 
problems and taking the limit but we do not include such details. We treat both bounded and unbounded 
domains together pointing out the differences when they appear. Let 

If the domain is bounded then by (j4} , 

^7 (t) = 2d J A{u)dx + 2 J J u{x)u{y)x ■ VJC(x - y)dxdy - J A(u)x ■ vdS 



= 2d J A{u)dx + J J (x-y)- V/C(x - y)u{x)u(y)dxdy - J A{u)x ■ v(x)dS, (54) 

where the second integral was obtained by symmetrizing in x and y, the time dependence was dropped for 
notational simplicity and v{x) denotes the outward unit normal of D at x G 3D. By translation invariance 
and convexity of D, we may assume without loss of generality that x ■ v{x) > 0. For the rest of the proof we 
may treat bounded domains and D — W 1 together, since for each, 

^Z(t) < 2d J A(u)dx + 2 J J u(x)u(y)x ■ VJC{x - y)dxdy. 

We use (B2) on K., to obtain 

j t I(t) < 2d J A(u)dx - 2d/pW{u) + CiM 2 . 

By (D3), (B4) and Lemma [ill 

A(u)dx = I A(u)dx + / A(u)dx 

J{u<R} J{u>R} 



< C{M) + (to - 1) / $(u)dx 

J{u>R} 

< C{M) + (to - 1) / $(u)dx. 



Using that 2d(m — 1) < 2d(m* — 1) = 2d/p we have, 

jl{t) < 2d{m - l)F{u) + C[M, Ci). 

We use the energy dissipation inequality ([7]) to bound the first term, 

^-I(t) < 2d(m - 1)J"(« ) + C(M, Cx). 
dt 

From this differential inequality, the second moment will be zero in finite time and thus the solution blows 
up in finite time if, 

C(M,d) 
2d{m — 1) 
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By Lemma we may always find initial data with any given mass M > such that this is true, since 
there exists infimizing sequences with vanishing second moments. The final assertion follows from Theorem 
SJ Indeed, we have 



7(0) 



2d(m - l).F(uo) + C(M, C ± ) 



□ 



6.2 Critical Case: Theorems [5] and [5] 

The proof of Theorem [5] follows the proof of Theorem [51 

Lemma 13. Define J^vi = {h £ i 1 n i°° : u > 0, = M\. Suppose tC satisfies (Bl) and A(u) satisfies 
(B3) for m > 1 and A > 0. Suppose further that the problem is critical, that is, m = m* and let M c 
satisfy (|9|). If M satisfies M > M c , then miy M T = — oo. Moreover, there exists an infimizing sequence with 
vanishing second moments which converges to the Dirac delta mass in the sense of measures. 

Proof. We may proceed as in the proof of Lemma [T2l but instead choose 9 € y(M c /M) 2 ~ m , l) . Let a = d/p. 
By optimality of C m * , as before there exists h* such that, 



0C m * < 



J J h*(x)h*(y)\x - y\ a dxdy 



\h*\\i~ m * \\h*\\™l 



<C„ 



(55) 



As above, we assume h* > and h* G 



Let n = \\h*\\\ /d M- 1 / d , A > and h\(x) = X d h* (Xfix) . By (Bl) and (B3), Ve > there exists a A and R 
sufficiently large such that by h* E C° 



"<OC 



\ dm — d ji/r 

(m* — l)||/r||i 

A a ii^ii?- m 1^n; 



-e)C m * / ll/i*!!! 



2 V M 

However in this case a = dm — d and m = m* , therefore by (I55[) and Lemma 21 



■F(/ia) < A 



dm —d\ \ l * 1 1 m 



h*\\l 



M(A + e) (0-e)C m * / " 2+a/d 



Then, 



-F(/ia) < A 



dm* —d 



(m* - l)\\h*\\! 2 

M(A+e) (6» — e) 



\h*\\%* 

\m\i 



(to* — 1) 



M 



C m *M 2 - a/d 



mi- 



Then since A/(m* - 1) = C m *M^- m /2 and a/d - 1 = 2 - to* we have, 



•F(/i A ) < A 



dm* —d 



\h*\ 



2m 



C m *M 2 - a/d 



A 



Mc 
M 



2-m 



Since > (M c /M) m we may take e sufficiently small and A — »■ oo to conclude that infy^ J 7 = — oo. As 
before, h\ converges to the Dirac delta mass in the sense of measures. □ 



Proof. (Theorem [5]) The theorem follows from a Virial identity as in Theorem [51 



□ 
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Proof. (Theorem [5]) As in Theorem [5] we have by (C2), (C3) and if D is bounded, the convexity of the 
domain, 

-^I(t) < 2dA J A(u)dx + J J u{x)u(y){x - y) ■ VK.(x - y)dxdy 
<2dAl(A-^+C 1 M 3 / 2 I 1 / 2 . 
Clearly, if M > M c then / — > in finite time if 1(0) is sufficiently small. □ 

7 Conclusion 

The prior treatments of (QJ have restricted attention to cither very singular kernels (PKS) or very smooth 
kernels (as in jB]). Moreover, most work has been restricted to power- law diffusion. We extend these 
approaches to develop a unified theory which applies to general nonlinear, degenerate diffusion and attractive 
kernels which are no more singular than the Newtonian potential. Existence arguments may apply to more 
singular kernels or unbounded initial data, however, to the authors' knowledge, Lemma [2] or something 
analogous must be available for any known uniqueness argument to hold. We generalize the existing notions 
of criticality for PKS and show that the critical mass phenomenon observed in PKS is a generic property of 
critical aggregation diffusion models. We extend the free energy methods of [23 [HI E] to obtain the sharp 
critical mass for a class of models with general nonlinear diffusion and inhomogeneous kernels. In particular, 
we show that the critical mass depends only on the singularity of the kernel at the origin and the growth of 
the entropy at infinity. The results presented here hold on bounded, convex domains for d > 2 and on R d for 
d > 3. 



8 Appendix 

8.1 Auxiliary Lemmas 

Lemma 14. Let F be a convex C 1 function and f — F' . Assume that f(u) G L 2 (0,T, H 1 (D)), u G 
H 1 ^, T, H^ 1 (D)) and F(u) G L°°(0,T, L l {D)). Then for almost all < s, r, < T the following holds: 

(F(u(x,t))-F(u(x,s))) dx = J («t, /(«(«))> dt. 

Lemma 15. Let F(u, t) G C 2 ([0, oo), [0, oo)) be a convex function such that F({)) = and F" > on (0, oo). 
Let f n , for n = 1,2,..., and f be a non-negative function on D bounded from above by M > 0. Furthermore, 
assume that f n — / inL 1 (D) and F(f n ) — > F(f) inL 1 (D), then |j/„ — fW^^ — > as n — > 0. 

Lemma 16 (Weak Lower-semicontinuity) . Let p e be non-negative L] oc (Dt) and f e a vector valued function 
in L} oc (D T ) such that V0 G C^°(D^)and^ G C* c °°(i5^, R d ) 

p e <j)dxdt — > I p<f)dxdt 

D T J D t 

f e ■ ^dxdt ->• / / • ^dxdt. 

D T J D t 



Then 



f - |/| 2 dxdt < liminf / — |/ e | 2 dxdt 
Jd t P Jd t Pt 
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8.2 Gagliardo-Nirenberg-Sobolev Inequality 

Gagliardo-Nirenberg-Sobolev inequalities are the main tool for obtaining LP estimates of PKS models and 
are used in many works, for instance [381 1131 1581 134] . The following inequality follows by interpolation and 
the classical Gagliardo-Nirenberg-Sobolev inequality. 

Lemma 17 (Inhomogeneous Gagliardo-Nirenberg-Sobolev). Letd > 2 and D C M. d satisfy the cone condition 
(see e.g. f^). Let f : D -> K satisfy f G L p n L q and V/ fe € L r . Moreover let 1 < p < rk < dk, 
k < q < rkd/(d — r) and 

Iks , 

< 0. (56) 

r q d 

Then there exists a constant Cqns which depends on s,p,q,r,d and the dimensions of the cone for which D 
satisfies the cone condition such that 

ll/lk. <C GJV 5||/DII/ fe ||^, (57) 

where < on satisfy 

1 = a\k + a 2l (58) 

and 

1 1 , — s lk, , 

= ai (— + . 59 

q p d r p 

Proof. We may assume that / is Schwartz then argue by density. Let /3 satisfy max(g, rk) < /3 < rkdj (d — r). 
First note by the Gagliardo-Nirenberg-Sobolev inequality, [Theorem 5.8, pQ], we have 

ll/ fc ||,3/fe <P,k,r,s \\f k \\l~ \\f k \\w°.r 

for /i G (0, 1) determined by interpolation and 9 = s~ 1 (d/r — dk/(3) G (0, 1). Moreover, the implicit constant 
does not depend directly on the size of the domain. Therefore, 

\\f% /k < i!/iir e)(i ^ )/(i - Mi - e)) ii/ fe ii^ Mi " e)) - 

Now, where A G (0,1) determined by interpolation, 

11/11, <||/||^ A) ||/ fc ||^ 

< ||^||(i-A)+(i-e)(i-p)/(i-/ 1 (i-e))||yfc||^/(fc-fcM(i-e))^ 

□ 



8.3 Admissible Kernels 

We now prove Lemmas 11121 and [3l We begin with the following characterizations of L p '°° . 

Lemma 18. Let Fix) = /(|a;|) G Lj oc n C° \ {0} be monotone in a neighborhood of the origin. If r~ d l p — 
o(/(r)) os r -4- 0, then F $ Lf£°. 

Proof. Since we have assumed / to be monotone in a neighborhood of the origin, without loss of generality 
we prove the assertions assuming / > on that neighborhood, since corresponding work may be done if / is 
negative. For any a > 0, by monotonicity, we have a unique r{a) such that f(r) > a,Vr < r(a). We thus 
have that A/ (a) = uidr(a) d , where uid is the volume of the unit sphere in R d . By the growth condition on / 
and continuity we also have that for a sufficiently large, 

K{a)- d/p < f(r(a)) = a. (60) 

Now, 

a p X f (a) =u d a p r{a) d . 
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□ 



Hence, by we have Ve > there is a neighborhood of infinity such that, 

io d a p r{a) d > e~* '. 

We take e ->■ to deduce that F ^ L p <°°. 

Lemma 19. Let F(x) = f(\x\) G i; 1 oc nC°\{0} 6e monotone in a neighborhood of the origin. Then f G Lfo^° 
?/ cmd only if f = 0{r~ d l p ) as r — » 0. 

Proof. Since we have assumed / to be monotone in a neighborhood of the origin, without loss of generality 
we prove the assertions assuming / > on that neighborhood. 

First assume that / ^ 0{r~ d l p ) as r — ► 0, which implies that for all 6q > and every C > there exists an 
r c < $0 such that 

f(r c ) > Cr c d/p . 

We now show that in a neighborhood of the origin, the function f(r) — Cr~ d l p is strictly positive for r < rc- 
Suppose not. Since both /, r~ d l p are monotone, there exists ro such that f(r) < Cr~ d / p for r < ro- However, 
this contradicts / =^ 0(r^ d ^) as r — > 0. Thus, we have that 

f(r) > Gr- d l p 

in a neighborhood of the origin (r < rc). Since for all C > we can find a corresponding rc, this is equivalent 
to r~ d ' p = o(/(r)), and by Lemma QJ] we have that / £ L p '°°. 

On the other hand, if / = <D{r~ d / p ) as r — >• there exists (5 > and C > such that for all r < S, 

f(r) < Cr' d/p . (61) 

By monotonicity, for all a > there is a unique r(a) G [0,(5] such that 

/(r) > a, for r < r(a), (62) 

where we take r(a) — if f(r) < a over the entire neighborhood. By (|6"Tj) and (|62l) . we have, necessarily 
that r{a) < oT p l d . Therefore, 

a p \ f (a) = a p cu d r(a) d < 1, 

which implies /l Sl(0 ) € L p '°°. □ 
Remark 10. Similar statements may be made about the decay of F(x) at infinity. 
Proof. (Lemma [lj By the fundamental theorem of calculus and condition (BD), 

/oo 
\d r d Xt d Xj JC(rx)\ dr 

<M- d . 

Similarly, this argument also implies |V/C| < |a;| 1 ~ d , which in turn implies V/C G x d /( d - 1 ). 00 . If d > 2 then 
we can carry out this argument another time and show that \IC\ < |a;| 2_<i . Moreover, in d = 2 we see that X] 
could have, at worst, logarithmic singularities at zero and infinity. 

□ 

Proof. (Lemma [2]) We compute second derivatives of the kernel JC in the sense of distributions. Let <f> G C^°, 
then by the dominated convergence theorem, 

/d x K.d x . (pdx — lim / d x K,d x 4>dx 
^°J\x\>e 

= - lim / d Xj K{x)^-Mx)dS - PV / d XzXj Kcj)dx. 

^ a J\ x \=e \ x \ J 
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By V/C G L d l l ~ d ~ 1 ^ oa and Lemma[H we have V/C = 0(\x\^ d ) as a; — > 0. Therefore for e sufficiently small, 
there exists C > such that, 



d Xj K,(x)f-Mx)dS 

x\=e \ x \ 



<Cl |ar| 1_,i |^(ar)|d5 



= C lex] 1 - 4 \<f)(ex)\ e^dS = C |0(O)| . 

J\x\ = l 

Similarly, we may define D 2 IC * <j> and we have, 

\\D 2 K * 0[[ p < CII^Ip + ||PV /" d XiX .K{y)<t>(x - y)di/|| p . 

Therefore, the first term can be extended to a bounded operator on LP for 1 < p < oo by density. The admis- 
sibility conditions (R),(BD) and (KN) are sufficient to apply the Calderon-Zygmund inequality [Theorem 
2.2 |56)]. which implies that the principal value integral in the second term is a bounded linear operator on 
L p for all f < p < oo. Moreover the proof provides an estimate of the operator norms, 



||PV / d XuXj K(y)u(x - y)dy\\ p 



</ p^tHIp 1<P<2 
p\\u\\ p 2<p<oo. 



□ 



Proof. (Lemma|3) The assertion that D 2 K £ implies /C £ -kf/ c 7 follows similarly as in Lemma[TJ 

Now we prove the reverse implication. Let K £ if/ c 7 2 ^°° w e s how that D 2 K = 0(r~ d / 7 ) as r —> 0. 
Assume for contradiction that D 2 /C ^ 0(r~ d li) as r ->■ 0. This implies that fc" ^ 0(r- d ^) or that fc'(r)r" 1 ^ 
0(r~ d /' 1 ) as r — >• 0. These two possibilities are essentially the same, so just assume that k" ^ 0(r~ d /' 1 ). By 
monotonicity arguments used in the proof of Lemma 1 1 91 this in turn implies r~ d ^ = o(k"). However, this 
means that for all e, there exists a 6(e) > such that for r £ (0, S(e)) we have, 



k{r) - k(S(e)) = - / k'(s)ds = / / k"(t)dtds + (r - 6(e))k' (6(e)) 

J 6(e) J 5(e) J 5(e) 

> e -\2-dh + ^ 

which contradicts the fact that k(r) = 0(r 2 ~ d < > ' 1 ) as r — > by Lemma [T9l 

The assertion regarding V/C is proved in the same fashion. □ 
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